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1. Introduction 

It is well known that there are many different criteria for the truth of the Riemann 
Hypothesis (RH). Some of these are not directly related to the important high level 
computations and developments concerning the non trivial zeros of the Riemann Zeta 
function. In fact, at the beginning of the century M. Riesz, and later G.H. Hardy and 
J.E. Littlewood (among other important results in number theory) found a criterion of 
"classical type" for the truth of the RH. The above criteria are related to some series 
involving values of the Zeta function outside the critical strip, i.e. at integers arguments 
of the Zeta function [7, 9], and in a numerical context, very accurate calculations are 
needed toward a "possible kind of verification" of the RH. 

In the literature important remarks have been given by leading mathematicians (see for 
example, those cited in [3]). We may think that such criteria may have a limited interest 
since, with them, one should work outside the critical strip. It is, in fact, true that in 
dealing with the above criteria one needs the use of arguments of the Zeta function 
outside the critical strip, and problems of interchange of summations are present. As an 
example, in the above criteria, if one uses the formula established by the authors, one 
should give a meaning to an integration over the real line, which exists only for finite 
intervals. In order to obtain finite numerical results which give "satisfactory" values to 
the functions supposed to be equal to the reciprocal of the Zeta function outside and 
inside the critical strip, the integration should be carried out using a special sequence 
of upper limit of integration extending to infinity [6]. 

But lately, there have been new developments and rigorous results in connection with 
this kind of problem: first a "regularization" of the representations of the Zeta function 
(a pioneering work by Maslanka [8]), followed (in particular) by a new rigorous discrete 
formulation with theorems concerning the above criteria (the works of Baez-Duarte [1, 

2, 3]). 

In light of these new approaches, we thought that some of the above criteria deserve still 
more study, at least in the direction of some numerical experiments. Thus, we introduce 
additional functions containing two parameters, in order to have additional confidence 
in the numerical results of the experiments. 

The content of this work is as follows: in Section 2 we define a general set of functions 
with two parameters a and /3 in the spirit of Riesz and of Hardy-Littlewood and then 
obtain the discrete "representation" of the Zeta function of our set by means of the 
two parameter Pochammer's polynomials with their coefiicients. For the reader the 
discussion of the conditions are then given in Appendix A and in Appendix B (they 
follow strictly the ingenious method of Baez-Duarte for the Riesz case a = P = 2). 
In Section 3 we then obtain in some "limit", a Poisson distribution for the coefficients of 
the Pochammer's polynomials; this is useful in the context of the numerical experiments. 
These are presented in Section 4 where many various limiting cases are treated. In the 
case of increasing values of the parameter f3, the experiments indicate that the Poisson 
distribution becomes more and more exact and the discrete function Ck becomes a 
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constant which can be evaluated. Finally, we present the results of an experiment 
carried out in the critical strip with different values of 9^(s). As a consequence, we may 
argue that in the context of the range of validity of the experiments we present, at large 
and at low values of k, the RH appears to be barely true. 



2. The model 



We now consider a set of functions with two parameters {a, (5) to obtain These 
are simply an extension of these two cases: the first (with a = /3 = 2) introduced and 
studied by Riesz [9], the second one (where a = 1 and /9 = 2) by Hardy-Littlewood [7]. 
Let /i(n) be the Mobius function of argument n, where 

{1, if n = 1 

(—1)'^, if n is a product of k distinct primes 
0, if n contains a square 



n=l 



Let So = p-\-it and s a complex variable. For 9^(s) > p — 1 one has ^ = ^ 
The two-parameters family of functions is given by: 

so that expanding the right-hand side in powers of x, we obtain: 



(1) 



-+1 



where 



1 ry-M!!)y-(zl)!£!,-(=?«)rf, 

-i)V 



dx 



k=0 



k\ C{a + Pk) 



(2) 



The function ip{x) was introduced by Riesz (case a = P = 2) and by Hardy-Littlewood 
(case a — 1, P = 2). 

If '0(x) ~ C.-P for some e and for large x, then 



oc — p _ 



\^{s)\ < 
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would exist and would eventually be given by ^ with C,{s) ^ if we choose 
Let p — \- For a — (5 — 2 we have: 
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On the other hand expanding (1) in a similar way, we have that: 
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r(-^) ^ '^11 I r / ' /3 

^ /? fe=0 r=l \ J 1^ 



Thus: 



V.(5) = 5]c,P,(^ + l) (3) 

fc=0 

where 11(1" -) are the Pochhammer polynomials and the functions: 

n=l ^ ' 

were already studied by Maslanka [8] and Baez-Duarte [2] in the special case a — (5 — 2. 
Let 9^(s) > p + e (e > and p e [1, oo[). From a theorem of Baez-Duarte [2], which says 
that I -Pit (s) I < A - A;~^(*) where A is a constant depending on for large values of k 
we have that: 

oo 

\^{s)\ <Y,\ck\k<'^^') (5) 

In Appendix A we show that if o; > 1 and /? > the following holds unconditionally: 

1 



k 



where 

1 / 1 ^ 

n=l ^ 



Then we obtain: 

\^{s)\ < ^-r < —A — < — ^ < oo 
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From this it follows that equation (6) below, which gives [C(s)] as: 

^ oo 

Vi^) = ^ = E + 1) (6) 

is valid for > 1. Now, still from the theorem of Baez-Duarte [2] i.e. that 



A 

< 



k ^ ^ 



and from the hypothesis discussed in Appendix B i.e. supposing the RH to be true for 
> p + e (e > and p G [i, oo[): 

I I B 

then the above series given by (3) converges uniformely. In fact for 9^(s) > p + e we 
have: 

°° 5 A °° C 

fc=0 '^'^ fc " fc=0 

Following Baez-Duarte the series ^^(.s) extends analytically ^ to the half plane 
9^(s) > p = |. We have thus obtained for our family of functions with parameters 
a, P that a necessary and sufficient condition for Ci^) ^ iii the half plane 9^(s) > p 
(p e [|,oo[) is given by: 

con s1" 

\ck{a,P)\< Ve>0,Va>l,V/3>0 (7) 

Remark 1 

Instead of using the Mobius function p in Cfc, one may use (for the numerical 
computations) the formula involving values of the Zeta function: 

Remark 2 

From the bound above it follows not only theoretically but also in the context of a 
numerical analysis that it will be equally difficult to treat the case p G [|, 1], for example 
p = |, as will be the case p = | + e with e small (see the last experiment concerning the 
critical strip where P = |, P = | and p = 1. Below in Section 4 we will also treat the 
case ot — \. 
Remark 3 

The condition for the truth of the RH using Riesz and Hardy-Littlewood functions ip{x) 
is essentially the same as the one using the discrete function Ck with k eN. In a previous 
work [6] independent of the present one (which essentially uses the Baez-Duarte idea 
and theorems) some numerical results were obtained for iIj{x) in the case of the Hardy- 
Littlewood function {a = 1, (3 = 2) by integration in the x-space. The discrete version 



Riesz criteria revisited 



6 



using the function of the discrete variable k [2, 8] has advantages in the numerical 
computations which will be presented below. Before this we present another way to 
control the function Ck in a numerical context. 

3. Poisson like distribution 

We still consider the function Cfc given by: 

k 

Ck 



n=l ^ ^ 



Then 



fJWkHi-^) 



n=l 
oo 



n=l 

^yl^^-J^^Aik,n,/3) 



n=l 

Notice that A < 0. For P large we set A = to obtain the following approximation: 

/x(n) 



oo 

Eu(n) __k 



Ck 

TV 

n=l 



In this approximation we see that Ck of (8) becomes equal to = /c) of (1). Moreover: 

r)^-fe 



n=l 

ii{n) ^ ( 1 ^ ^ 



n=l p=0 ^ ^ 



Thus: 



p=Q ^' 

We are in the presence of a Poisson distribution for the c^: in this way, in our numerical 
computations, we may control in a "more satisfactory" way the values of c^. The 
approximation for Ck by means of the Poisson distribution for the c^'s we found, will be 
more satisfactory with increasing values of and for large values of k. We will also use 
the approximation given by (9) in which the upper limit of summation will be given by 
2k instead of oo, i.e. for large k, 

^^ = E^^^^"' (10) 

p=Q ^' 



Riesz criteria revisited 



7 



4. Numerical experiments 
4.1. The case a = | and (3 = 4: 

This is a case of interest since the behaviour of the Ck at large values of k is expected 
to be the same as the case a = /? = 2 [2, 8]. In fact from (7) we ask that for > |: 

|c,(7/2,4)| < ,,,,,,,,, ~ ~ |cfc(2,2)| (11) 
k 4 ki 

In the figures 1, 2 and 3 we give the plot respectively of log |cfc|, log(|cfc log fc|) and 
log(|cA;(log fc)^|) as a function of logfc for k up to 1000 together with the straight line 
with slope — I which is tangent to the curves at some point. The Ck were computed 
calculating (4) until n = 10000. 

This experiment indicates that Ck decays more fast than 3 as announced by Baez- 

(logfe)^fc3 

Duarte in [3] for the case a = P = 2, i.e. more fast then the bound (11) if the RH is 
true (see the necessary condition in Appendix B). 



log ] Cjt [ 

-2- 




log ( ] Cjt [ log k) 




Figure 2. Plot of log(|cfc| logfc) = C — | logfc together with the tangent straight Hne 
of slope — I . 
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log(|c^[ (loqk)') 




log k 



Figure 3. Plot of log(|cfe|(logA;)'^) = C — |logfc together with the tangent straight 
Une of slope — |. 



4.2. The case a = ^ 

If a = then all the Cfc(a, /?) are expected to have a decay similar to (7). Below 
(see figure 4) we present the numerical results for (3 from 1 to 6. This is in agreement 
with the theory of Appendix B. 
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Figure 4. Plot of Ck(^^^,P) for (3 — 1,2,3,4,5,6 (dotted lines) vs. the reference 
curve /(fc) — — 0.4fc~4. 



4.3. The case a = | with [3 ^ 00 

Let a = I be fixed, from (8) as [3 increases we get: 



i (k\ 1 



^lim^c, =;im^E(-iyO)c(7/2+,,) 



\0) C(7/2) 



i + i + E(-)(-ir 



feT-i + Ee)(-i: 

i=o 



C(7/2) 



Thus: 



1 



lim cu = , , X 
C(7/2) 



-0.112479 



VA; G N 



(12) 
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Figure 5. Plot of Ck for a = 7/2 and /? = 4, 5, 6, 7, 20 (from top to bottom). 



Our numerical experiments convalidate these results. We calculated the first 100 Ck for 
/? = 4, 5, 6, 7, 20. For /3 = 20 we get already a convergence to the theoretical limit (12), 
see figure 5. So, this infinite /3 limit obtained by the numerical calculations for low 
values of k (up to 100) indicates that RH may barely be true (see (7) as ^ —>■ oo). 



4.4- The Poisson distribution 

To demonstrate the goodness of the approximation's formula (10) we computed the Ck 
until k — 1000 for the case a — |, /3 = 4 (using (4)). Then using these already computed 
Ck we calculated also the first 500 Ck of (10). We ploted these two curves together. In 
figure 6 we see that from A; = 40 the Poisson approximation is essentially the same as 
the real function. 



4-5. A case with different values of p in the critical strip 

This experiment is carried out in the critical strip in order to compare the behaviour of 
Ck in three different cases of p. In fact from Appendix B the Ck should all decay at least 
as: 

\Ck\ < -T— (13) 

In particular for the cases we treat, i.e. a = /5 = 2 (the case of Riesz at p = |), 
a — (5 — 'iaXp = \ and a = P — Aatp — 1, the results are presented on figure 7, 8 
and 9 respectively. The plots indicate similar behaviour not in disagreement with (13). 
Further, by looking at these plots, the numerical computations visualize that: 

|CA:(2,2)| < |cfc(3,3)| < |cfc(4,4)| 

We know that |cfc(4,4)| < A;"i+^ for p > 1 and so |cfc(3,3)| -C k~i^\ This would 
indicate that there is no zero for p > |. Finally |cfc(2,2)| <^ k~4+^ would confirm that 
there is also no zero for p > |. 
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Figure 6. Plot of Ck for a = = 4 (left) vs. the Poisson approximation (right). 
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Figure 7. Plot of log(|cfe|) for a — P — 2 together with the straight line of slope 
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log 1 Ck [ 




4.6. The case a = | and [3 = 1 

In this case, from the Riesz formulation (1) or from (6), [C('5)]~^ is different from infinity 
if Ck decays at least as . Setting x = k the Riesz function (1) is given by: 

n=l 

Below (figure 10) we present the plot of /(fc) together with that of g{k) = — 0.07/c~^ for 
k up to 600 and for a maximum value of only some hundred for n. 

4. 7. The case a = \ 

In this case (a < 1!) we cannot employ the argument of Appendix A, but we have for 
^{s)-t>p= \: 

00 00 
|^(.)| < 5^A:-(^+^) |c,(l/2,/3)| < |c,(l/2,oo)| 

fc=0 fc=0 

From Subsection 4.3 we know that |cfc(l/2, oo)| = | ^(3/2) — 1| — 0.617, thus for any finite 
/?, V9(s) is also finite. 
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f {left) and g [right) 
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Figure 10. Plot of / and (a = |, /? = 1). 



5. Conclusions 

In this work we have revisited the original criteria of Riesz and of Hardy and Littlewood 
for the Riemann Hypothesis in light of recent pioneering works concerning the possible 
representations of the Riemann Zeta function by means of the Pochammer's polynomials. 
The discrete representations in the case a = [3 = 2 are do to Maslanka and to Baez- 
Duarte. In order to carry out our numerical experiments related to the criteria, we have 
first extended the analytical formulation to a more general class of functions containing 
two parameters a and /?; using a theorem of Baez-Duarte we have specified a sufficient 
and necessary condition for the truth of the RH for our general class of functions i.e. 
for the decay of the coefficients Ck as a power law of k. Moreover in doing this we have 
found the emergence of a Poisson-like distribution for the which should be exact in 
the large [3 limit. Numerical experiments have been carried out for various cases. 

(i) For a = I and /3 = 4 we have presented intensive calculation using the Mobius 
function up to = 10000 and for k up to some hundreds. For this case, the 
power law decay is the same as that appearing in the original work of Riesz 
(a = /? = 2) and also investigated numerically by Baez-Duarte. The experiments 
confirm the correctness of the power law within the range of the values of n and of 
k we were able to treat. 

3 

(ii) For a and j3 such that the Ck should all give the power law decay k~^ at large values 
of k to ensure the truth of the RH, i.e those where have presented some 
experiments at low values of k and some values of /3 which confirms this power law 
decay. All functions assume negative values with plots lying above a fixed curve 

3 

of equation y = Ak~i for some fixed constant A independent of /3. In addition for 
this expected behaviour we have reported the results of intensive computations for 
the cases a = (3 = 2,a = P = 3 and a = /? = 4. 

(iii) Finally we have considered some experiments in the large /3 limit which indicate 
that the plots of Ck become more and more fiat, well approximated by the mean 
value of the Poisson-type distribution we have founded. As (3 becomes large and 
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large the approaches a constant value, for all k, indicating that in this sense the 
RH may barely be true. 

This work, still accompanied by numerical experiments, may be expanded in the 
search of other new representations of the Riemann Zeta function, different of the one 
considered here; moreover there is the aim that the new criteria will be useful in the 
context of more numerical experiments. The works will be presented in a near future 
[4, 5]. 



Appendix A 

We follow strictly the lines of calculations of Baez-Duarte [2] to show that the 
representation (6) for [Cl^)]"^ is unconditionally valid for lH(s) > p — l,a > 1 and 
P > 0. We consider the quantity: 

n=l ^ 

Using the MacLaurin series (restricting ourselves to the main contribution), we have 
that: 

Then with the variable change y = ^ we obtain: 



1 k+l-l 

gj^y ' dy 



^k ^ 

-B{^,k + l) 

where 

Bix,^) = J^x (i-x) ^-=r(XT^ 

is the Beta function. 
Thus for k large, we have 

^ 1 ,q; — 1, A; ^ 1 

P P k ^ k 1^ 



Appendix B 

Still following Baez-Duarte [2] and here for the family of functions with parameters a 
and (3, we now show the necessity of the condition (7), assuming the RH to be true in 
the seminfinite strip IH(.s) > p — ^. 
We set M{x) — ^ P'{ri), then we obtain Ve > 0: 

n<x 



M{x) < X' 
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Integration by parts gives for the main contribution: 

With the variable change y = ^, using M{-^) <C y"^"" for y i (RH) we have: 

Jo Jo 
and finally with y^ = z we obtain 

a a-p-e , ^ fc+1 — 1 a-p-e+13 -. 

\ck\ < - z 1^ dz + k z ^ (l--^) dz 

P Jo Jo 

which for large k is given by: 

\Ck\ < '■ 1 -^3 < 

P k ^ k & ^ k ^ 
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